Concussion Prone Scenarios: A Multi-Dimensional Exploration in Impact Directions, Brain Morphology, and Network Architectures Using Computational Models.

While individual susceptibility to traumatic brain injury (TBI) has been speculated, past work does not provide an analysis considering how physical features of an individual's brain (e.g., brain size, shape), impact direction, and brain network features can holistically contribute to the risk of suffering a TBI from an impact. This work investigated each of these features simultaneously using computational modeling and analyses of simulated functional connectivity. Unlike the past studies that assess the severity of TBI based on the quantification of brain tissue damage (e.g., principal strain), we approached the brain as a complex network in which neuronal oscillations orchestrate to produce normal brain function (estimated by functional connectivity) and, to this end, both the anatomical damage location and its topological characteristics within the brain network contribute to the severity of brain function disruption and injury. To represent the variations in the population, we analyzed a publicly available database of brain imaging data and selected five distinct network architectures, seven different brain sizes, and three uniaxial head rotational conditions to study the consequences of 74 virtual impact scenarios. Results show impact direction produces the most significant change in connections across brain areas (structural connectome) and the functional coupling of activity across these brain areas (functional connectivity). Axial rotations were more injurious than those with sagittal and coronal rotations when the head kinematics were the same for each condition. When the impact direction was held constant, brain network architecture showed a significantly different vulnerability across axial and sagittal, but not coronal rotations. As expected, brain size significantly affected the expected change in structural and functional connectivity after impact. Together, these results provided groupings of predicted vulnerability to impact-a subgroup of male brain architectures exposed to axial impacts were most vulnerable, while a subgroup of female brain architectures was the most tolerant to the sagittal impacts studied. These findings lay essential groundwork for subject-specific analyses of concussion and provide invaluable guidance for designing personalized protection equipment.

Brain architecture-based vulnerability to traumatic injury.

The white matter tracts forming the intricate wiring of the brain are subject-specific; this heterogeneity can complicate studies of brain function and disease. Here we collapse tractography data from the Human Connectome Project (HCP) into structural connectivity (SC) matrices and identify groups of similarly wired brains from both sexes. To characterize the significance of these architectural groupings, we examined how similarly wired brains led to distinct groupings of neural activity dynamics estimated with Kuramoto oscillator models (KMs). We then lesioned our networks to simulate traumatic brain injury (TBI) and finally we tested whether these distinct architecture groups' dynamics exhibited differing responses to simulated TBI. At each of these levels we found that brain structure, simulated dynamics, and injury susceptibility were all related to brain grouping. We found four primary brain architecture groupings (two male and two female), with similar architectures appearing across both sexes. Among these groupings of brain structure, two architecture types were significantly more vulnerable than the remaining two architecture types to lesions. These groups suggest that mesoscale brain architecture types exist, and these architectural differences may contribute to differential risks to TBI and clinical outcomes across the population.

An interdisciplinary computational model for predicting traumatic brain injury: Linking biomechanics and functional neural networks.

The brain is a complex network consisting of neuron cell bodies in the gray matter and their axonal projections, forming the white matter tracts. These neurons are supported by an equally complex vascular network as well as glial cells. Traumatic brain injury (TBI) can lead to the disruption of the structural and functional brain networks due to disruption of both neuronal cell bodies in the gray matter as well as their projections and supporting cells. To explore how an impact can alter the function of brain networks, we integrated a finite element (FE) brain mechanics model with linked models of brain dynamics (Kuramoto oscillator) and vascular perfusion (Balloon-Windkessel) in this study. We used empirical resting-state functional magnetic resonance imaging (MRI) data to optimize the fit of our brain dynamics and perfusion models to clinical data. Results from the FE model were used to mimic injury in these optimized brain dynamics models: injury to the nodes (gray matter) led to a decrease in the nodal oscillation frequency, while damage to the edges (axonal connections/white matter) progressively decreased coupling among connected nodes. A total of 53 cases, including 33 non-injurious and 20 concussive head impacts experienced by professional American football players were simulated using this integrated model. We examined the correlation of injury outcomes with global measures of structural connectivity, neural dynamics, and functional connectivity of the brain networks when using different lesion methods. Results show that injurious head impacts cause significant alterations in global network topology regardless of lesion methods. Changes between the disrupted and healthy functional connectivity (measured by Pearson correlation) consistently correlated well with injury outcomes (AUC≥0.75), although the predictive performance is not significantly different (p>0.05) to that of traditional kinematic measures (angular acceleration). Intriguingly, our lesion model for gray matter damage predicted increases in global efficiency and clustering coefficient with increases in injury risk, while disrupting axonal connections led to lower network efficiency and clustering. When both injury mechanisms were combined into a single injury prediction model, the injury prediction performance depended on the thresholds used to determine neurodegeneration and mechanical tolerance for axonal injury. Together, these results point towards complex effects of mechanical trauma to the brain and provide a new framework for understanding brain injury at a causal mechanistic level and developing more effective diagnostic methods and therapeutic interventions.

A Multibody Model for Predicting Spatial Distribution of Human Brain Deformation Following Impact Loading.

With an increasing focus on long-term consequences of concussive brain injuries, there is a new emphasis on developing tools that can accurately predict the mechanical response of the brain to impact loading. Although finite element models (FEM) estimate the brain response under dynamic loading, these models are not capable of delivering rapid (∼seconds) estimates of the brain's mechanical response. In this study, we develop a multibody spring-mass-damper model that estimates the regional motion of the brain to rotational accelerations delivered either about one anatomic axis or across three orthogonal axes simultaneously. In total, we estimated the deformation across 120 locations within a 50th percentile human brain. We found the multibody model (MBM) correlated, but did not precisely predict, the computed finite element response (average relative error: 18.4 ± 13.1%). We used machine learning (ML) to combine the prediction from the MBM and the loading kinematics (peak rotational acceleration, peak rotational velocity) and significantly reduced the discrepancy between the MBM and FEM (average relative error: 9.8 ± 7.7%). Using an independent sports injury testing set, we found the hybrid ML model also correlated well with predictions from a FEM (average relative error: 16.4 ± 10.2%). Finally, we used this hybrid MBM-ML approach to predict strains appearing in different locations throughout the brain, with average relative error estimates ranging from 8.6% to 25.2% for complex, multi-axial acceleration loading. Together, these results show a rapid and reasonably accurate method for predicting the mechanical response of the brain for single and multiplanar inputs, and provide a new tool for quickly assessing the consequences of impact loading throughout the brain.